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Abstract 

In order to clarify the puzzling problems in double charm production, relativistic corrections at 
order v"^ to the processes e^e~ — > J/V' + Vc and e"^e~ — > J/il^ + cc at B factories are studied in 
non-relativistic quantum chromodynamics. The short-distance parts of production cross sections 
are calculated perturbatively, while the long-distance matrix elements are estimated from J/^lJ 
and 7]c decays up to errors of order v^. Our results show that the relativistic correction to the 
exclusive process e^e~ J/ij) + r]c is significant, which, when combined together with the next-to- 
leading order as corrections, could resolve the large discrepancy between theory and experiment; 
whereas for the inclusive process e^e~ J /ip + cc the relativistic correction is tiny and negligible. 
The physical reason for the above difference between exclusive and inclusive processes largely lies 
in the fact that in the exclusive process the relative momentum between quarks in charmonium 
substantially reduces the virtuality of the gluon that converts into a charm quark pair, but this 
is not the case for the inclusive process, in which the charm quark fragmentation c — > J/ij: + c \s 
significant, and QCD radiative corrections can be more essential. 
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I. INTRODUCTION 



The exclusive and inclusive double charm production in e"'"e~ annihilation at B factories 
has been one of the most puzzling problems in nonrelativistic QCD (NRQCD) [l| and heavy 
quarkonium physics for years. The measured cross sections for exclusive double charmo- 
nium production e~^e~ J/ip + Vc and inclusive double charm production e~^e~ J/ip + cc 
at -^i = 10.6 GeV are much higher than the leading order (LO) predictions in NRQCD. 

The cross section of e+e~ ^ J / if) + rjc aX y/s = 10.6 GeV observed by Bellejs, 4] is 



(7[e+e- ^ J/^ + r]c] x 13'''[> 2] = 25.6 ± 2.8 ± 3.4 fb, 



(la) 



and by BaBar[5] is 



cT[e+e~ ^ J/ip + ric] x i3''=[> 2] = 17.6 ± 2.8 ± 2.1 fb, 



(lb) 



where B''^''[> 2] is the branching fraction for rjc decay into at least two charged particles. 
Since B'^''\> 2] < 1, Eq.[T] is the lower bound of the cross section of the double charmonium 
production. It is about an order of magnitude larger than the theoretical predictions^, Q, Q], 
at the leading order of strong coupling constant and quark relative velocity v in the 
NRQCD factorization approach|l|. Injp] the next-to leading order (NLO) QCD result is 
given and it is found that the 0{as) corrections may greatly reduce the large discrepancy 
between theory and experiment. Meanwhile, the relativistic corrections are also considered 



by a number of authors. In 



lOj, relativistic and bound state effects are discussed on the 



basis of relativistic quark models. In 
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15| . relativistic effects are estimated 



in the framework of the light cone method, and very large enhancement effects on the cross 
sections can be found with certain light-cone distribution amplitudes. In [l^, the authors 
calculate the relativistic correction based on NRQCD and light cone method with long- 
distance matrix elements determined from the Cornell potential model, and they point out, 
however, that after subtracting parts of the light-cone distribution functions that correspond 
to corrections of relative-order in NRQCD, the enhancement effect due to relativistic 
corrections is not very large but still substantial. Probably, with both relative-order-a^ and 
relativistic corrections we may resolve the discrepancy between theory and experiment for 
e+e- ^ J/^ + r/eQ . 

The large discrepancies between theory and experiment exist not only in the exclusive 
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double charmonium production processes such as e'^e — > J/ip + r]c but also in the inclusive 
production process of e^e^ ^ J/tp + cc. The experimental result measured by Bellelsj] is 

a[e+e- ^ J/ip + cc + X] = O.SrtoH ± 0.17pb, (2) 

which is about five times larger than the theoretical calculations [itI. [isl . [lol . 20, 21| based 
on NRQCD at leading order of a,^ and v'^. In [22] the two-photon contribution and in 2l| 
the color-octet contribution to J /ip + cc production are further considered, but they are too 
small to resolve the large discrepancy. Other attempts to solve the problem can be found in 
a comprehensive review on heavy quarkonium physics [2]. It is certainly interesting to see 
whether the large discrepancy can be resolved by inclusion of NLO QCD corrections and 



relativistic corrections. Recently, in j23|, the authors calculate the prompt J/i}) + cc + X 
production at NLO including direct production and feeddown contribution mainly from 
the il){2S), and find the NLO as corrections are large and positive, and could be helpful to 
settle the problem between experiment and theory. But we still need to know how large are 
the relativistic corrections, and whether they are positive or negative to the solution of the 
problem. 

In this paper, we consider the relativistic corrections to both these exclusive and inclusive 
double-charm production processes based on NRQCD formulas in the color-singlet sector, 
since the color-octet contributions are negligible |2l|. In order to avoid ^ the model depen- 
dence in determining the long distance matrix elements, differing from [16], we determine 
the matrix elements of up to dimension-8 four fermion operators from the observed decay 
rates of J/ip and rjc- We find that the relativistic effect on the double charmonium pro- 



duction e^e'^ —>■ J/ip + rjc is substantial and comparable to the estimate of 16|; whereas 
for the inclusive production e^e^ —>■ J/ip + cc the relativistic corrections are very small 
and negligible. The rest of this paper is organized as follows. In Section II, we give the 
general formulas of the production rates in NRQCD at order. Relativistic corrections 
to the exclusive process and the inclusive process are studied in Section III and Section IV 
respectively. A summary for the f ^ order corrections to these two processes will be given in 
Section V. In the Appendix we give some of the analytic results. 
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II. PRODUCTION CROSS SECTIONS IN NRQCD 



In NRQCD the production and decay of charmonia are factorized into two parts, the 
short distance part that can be calculated perturbatively, and the long distance part can 
be estimated by lattice calculation, phenomenological models, or from other experimental 
observables. The long distance parts are related to the four fermion operators, characterized 
by the velocity v of the charm quark in the meson rest frame. The production cross sections 
of rjc and J/i}) up to order are jl] 



■mt 



nil 



FA'S, 



mt 



m'Stm + ^i^^(0|nesf )|0) + 0{v'a) 



The operators are defined as 



2 ' ^ " '2 

The short distance coefficients can be evaluated by the matching condition: 



pert NRQCD 



(3a) 
(3b) 

(4a) 

(4b) 
(4c) 
(4d) 

(5) 



The left hand side of Eq. 5] can be calculated by the spinor projection method 24 1. 



Furthermore (see, e.g., 



251] [26|]), the projection of v{P/2 — q)u{P/2 + q) onto a particular 



angular momentum state can be expressed in a Lorentz covariant form. The momenta of 

n 

quark and antiquark in an arbitrary frame are respectively [27|: 



ip + g = L(^P, + q) 



(6a) 



-P -q = Li-Pr 
2 ^ 4 



q), (6b) 

where P^ = {2Eq, 0), Eg = ^/m^~+~^, and 2q is the relative momentum between two quarks 
in the meson rest frame. is the boost tensor from the meson rest frame to an arbitrary 
frame. 
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In the meson rest frame the expression of projection onto a state of 5" = is 

v{-q, X2)u{q, Ai)(i, Ai, A2IO, 0) = 

A1A2 

In an arbitrary frame it becomes 

Y v{q, X2)u{q, Ai)(2> Ai, -, A2IO, 0) = 

A1A2 

11 -\- IE 1 

' /'-^-m)75^^-^(-/' + ^ + m). (8) 



2y/2{E^myr ^ '"^ 2E '2 
Similarly, expressions of projection onto a state of = 1 in the rest frame of the meson 
and an arbitrary frame are: 

v{-q, AaMq, Ai)(^, Ai; ^, A2II, e) = 



A1A2 



Y v{q, X2)u{q, Ai)(^, Ai; ^, Asjl, e) = 
A1A2 

In our calculation the Dirac spinors are normalized as uu — —vv — 2mc. Then the short 
distance part of the cross section can be calculated at any order of v. 

At order v"^ the cross section of e^e' J/ip + ric and e~^e~ — >• J/ip + cc can be expressed 

as 



(0|Pi(35f)|0) ,r,irn (lcVc^\r,\ 1 



3 

(0|yi('S;-)|0> ^ jMi,. d L/PS . (lOa) 



/7^dL/P5+<°l^'';^*'l«>l/iV.<iL/P5. (10b) 



cr(e'''e J/ip + cc) 

(o|0i(3^f)|o) 1 , {ow.esfm 1 

3 2s 

Here LIPS means the Lorentz invariant phase space, bar means averaging spins over the 
initial states and summing spins over the final states, and Mq, Mu/^, Mi^^, Nq, Ni can be 
calculated perturbatively. 



+2 flipped diagrams 



FIG. 1: e+e- ^ {00)33, + (cc)i5o 
III. RELATIVISTIC CORRECTIONS TO e+e ^ J/ip + r]^ 
A. Short distance part 

There are four Feynman diagrams in the process e+e~ —>■ {cc)ss^ + [cc)iSo, shown in Fig. [1]. 
The amphtude of the process can be expanded in terms of the quark relative momentum in 
charmonium: 

M(e+e- ^ (cc)3s, + (cc)iso) = {'^'^y^^A{q^,qJ = 



+ 



1 , ^ d'A 



Qii—Qvc 



El E2 

^ a d^A 



=0 2 



dq^dqi 



<ii:=qvc=o 



and ^(g',/., (?r)J is expressed as 
A{q^,q^J = 



- Yl I](^,Ai;^,A2|l,5,)(^,A3;^,A4|0,0)(3,^;3,j|l)(3,fc;3,/|l) 

A1A2A3A4 ijkl 

PI PI P2 P2 



(11) 



(12) 



where (3,z;3,j|l) = 5ij/\fWc and (3, A;;3,/|l) = 5ki/\fN'c are the color-SU(3), Clebsch- 
Gordon coefficients for QQ pairs projecting onto a color singlet state. Using Eq.[8] and 
Eq. [9b], we can express A{q^,qr^^) in a covariant form. The factor 



^c'^^Cn 1/2 r 

) ' comes from 

El E2 



the relativistic normalization of the cc state, and Ei = ^yrri^ + q^^^ jr;^ = ^/m^ + q,,^^. 
For the 5* wave charmonium q'^q^ = |q^(— fif"^ + ^-^-), where P^ = AE'^, P ■ q = and 
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the odd-power terms of v vanish. Then at leading order of f ^ 



^?;?J^r^^(0, 0) + l:qlqlA;:^A{0, 0), 



(13) 



where A^^ = A^^ = 

To obtain the short distance part E = + should also be expanded in f ^ = 

Then the relation of Mq, Mu/^, Mi^^ to |Mp can be easily written down. Details of how to 
get the short distance coefficients from the covariant projection method can be found in e.g. 



28| 



Mn 



:{A{0,0)A*{0,0)) 



(14a) 



1 



a(^A(0,0)A*(0,0)) 



<=0 



+^{a;^u^p{Pj/^)a*{o, 0) + a;:'^u^p{Pj/^)a{o, o))) 



(14b) 



M 



a(£-A(0,0)A*(0,0)) 



(2Ar,m,) 



1=0 



+-(A;;f n„^(p,jA*(o, 0) + n„^(p,jA(o, 0))) 



(14c) 



Here Mq is just the leading order result, and it agrees with the previous result in ref.j?] 



_ 1 2048(s- 16m2)(l + cos2(^))(47ra)2(47ra,)2e 



2^2 
1 



9m 



(15) 



The expressions of Mu/^ and Mi^^ are 



1 + cos^) {AnayiAiTasrel 512{lQml - 3s) 512(2560m^ - 592mls + 27s 



+ 



(16a) 



(l + cos^)(47ra)2(47ra,)2e^ _16384 1024(lls - 80m2)(s - 16m 



M 



l»7c 



9m 



+ 



3m^s^ 



(16b) 



where Nc = 3, = |, and 6* is the angle between the J/ip and the electron. 
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B. Long distance part 



In this section we present a phenomenological estimate of the color-singlet production 
matrix elements which are extracted from charmonium decay data. It is known that up to 
errors of order the color-singlet production matrix elements are related to the decay matrix 
elements through vacuum saturation [1]. There are three independent NRQCD matrix 
elements at order v^, (0|Ci(i5o^^)|0), (0|Oi(35f )|0), and (0|Pi(i5o^^)|0) = (0|Pi(35f )|0)/3 
(1 + 0(f^)). We can get them through the J /iIj leptonic decay J /ip ^ e~^e~ and hadronic 
decay J/ip ^ LH, and the rjc photonic decay rj^ — > 77. The theoretical results at next-to- 
leading order of and for J/ip ^ e'^e~, J/ip ^ LH ^, and t]c 77 are summarized in 



ref.|25| as 

n,. = 2eJW((l - (2°-^>.)(°|O.Cf)l°> _ 1 (W^fMy 

V[J/^ e+e-] = _ 16c^.^ (O|0i(^gp|O)/3 _ 4 (0|n(^^p|0)/3 x ^^^^^ 

243 V TT 127r^ — 108 / 

(17c) 

Solving these equations at leading order of as (QCD radiative corrections not included), we 

get 

(0|Oi(^S'o^^)|0) = 0.286G'ey^ (0|Oi(3^f )|0)/3 = 0.295Gey^ (18a) 
mlpR^m^^,.,n.l,-^GeV\ (18b) 
for rrir =1.5Gev and = 0.26. The experimental data of these decay rates can be found 



from [29|], and we choose their central values T[J/ip e+e^] = 5.55KeV, T[J/iIj — > LH] 
69.3KeV, and T[ric —>■ 77] = 7.14KeV. The matrix elements can be expressed as functions 
of the charm quark mass, 

(0|Oi(iS;'^)|0) = 0.127m2, (19a) 
mMM^O.m,nl (19b) 



We do not include the electromagnetic process oi J/ip ^ -f* LH. 
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<»l^-''f»l°> = <»l^-'t'l°>=0.014.ng (19c) 

Including the QCD NLO radiative corrections, and doing the calculation in the same way 
as above, we have 

(0|C»i(^^o^^)|0) = 0.432Gey^ (0|Oi(3^f )|0)/3 = 0.573Gel^^ (20a) 

0.514 X 10"^G'e1/^ (20b) 



(0|Pi(^gg-)|0) _ (0|Pi(^gf)|0) 
ml 3m^ 



for rric =1.5Gev and Us = 0.26. And we then have 

{0\OiCS^')\0) = 0.192ml (21a) 

<^MM = 0.255»^, (21b) 

3 

<^^iM=<»!>>= 0.023< (210 

as functions of m^. 

n 

quark pole mass and the quarkonium mass, and in pj| the matrix elements are calculated 
based on the potential model. Differing from their methods, we get the production matrix 
elements by using the experimentally observed charmonium decay rates. As shown in[l|, 
the difference between the color-singlet production and decay matrix elements are of order 
v^. So our method should be valid at order f^, and our estimates of the production matrix 
elements should be good numerically, if the high order QCD and f ^ corrections for the decays 
are small and the experimental errors are not large. 



C. Numerical results and discussions 



At ^/s = 10.6GeV, with mc = l.SGeV and as = 0.26, using the matrix elements in 
Eqs. |18lfT9] determined from charmonium decays with v"^ corrections but without QCD ra- 
diative corrections, and making the phase space integral, the leading order cross section of 
e'^e~ J/ip + rjc (LO means for the short-distance part, since for the long-distance ma- 
trix elements the corrections are already included) is 3.07 fb. The relativistic correction 
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f 3.0- 



2.5-. ' " ' 

2.0- 

1.5- 

1 1 1 1 1 

1.4 1.5 1.6 

m^(Gev) 

FIG. 2: e^e~ J/ip + rjc cross sections with relativistic corrections to long-distance matrix 
elements extracted from charmonium decays (without NLO QCD radiative corrections). The lower 
line represents the LO result in v, and the upper line represents the result with v'^ corrections to 
the short-distance coefficients. Here the coupling constant is fixed as as = 0.26. 

contributes 0.798 fb, which gives about 26% enhancement and results in 3.87 fb for the 
cross section. The cross sections as functions of rric are shown in Fig[2]. The lower line 
represents the LO result in v, and the upper line represents the result with f ^ corrections. 
Since the long-distance matrix elements are proportional to the squared quark mass, and 
the short distance coefficients are found to be quite stable when rric changes, so the cross 
section goes down as becomes smaller. However, these results are obtained for fixed 
value of as = 0.26, and if the running of = as{2mc) is assumed, the rric dependence of 
cross sections will be substantially changed, making the cross section at mc=lA GeV close 
to that at mc=1.6 GeV. 

If using the matrix elements in Eqs. [20|2T] determined from charmonium decays with 
both and as corrections, at order f° the cross section is 9.02fb, and at order v"^ the cross 
section is 11.26fb for as = 0.26 and rric = l.SGev. When rric varies from 1.4Gev to 1.6Gev, 
the cross section goes up from 9.18fb to 13.43fb. 

The experiment result in Eq. [1] is about an order of magnitude larger than the leading 
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1.4 



1.5 



1.6 



m(GeV) 



FIG. 3: e+e^ J/ip + rjc cross sections with relativistic corrections to long-distance matrix 
elements extracted from charmonium decays (with NLO QCD radiative corrections). The lower 
line represents the LO result in v, and the upper line represents the result with v'^ corrections to the 
short-distance coefficients. Here the coupling constant is fixed as Og = 0.26. Note that the QCD 
radiative corrections to the short-distance coefficients (with K=1.8) are included for the upper line 
but not the lower line. 

order result . The result at order f ^ shows that relativistic corrections can enhance 

both the short distance coefficients and the long distance matrix elements. The next-to- 
leading order v"^ coefficient (Mu/^ + Mi^^)m^ is about 2.2 times larger than M for = 
l.SGev and a, = 0.26, and the matrix elements of {0\Oi(^Sli')\0) and {0\Oi{^Sf)\0) also 
become about 1.17 times larger than that in the leading order calculation. These make the 
cross section become 1.7 times larger after including the relativistic effect. If we determine 
the matrix elements including the QCD radiative correction, at f ^ order the theoretical 
result is 11.26fb for as = 0.26 and rric = 1.5Gev. It shows that the relativistic corrections 
are significant. 

In the above discussions, we have not considered the QCD radiative corrections to the 
short- distance coefficients as shown in tefl. If we further combine the NLO QCD corrections 
(9| with the corrections with a fixed value = 0.26, then the cross section will go from 
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TABLE I: Experimental and calculated cross sections of a[e~^e J/ip + rjd with nic =1.5 GeV 
and as=0.26. See text for the definitions oi aLo(as,v2), (^NLO(as)^ (^NLOiv^)^ and (Tj^Loias,v2)- 



Experimental Result 



J/^ + ilc] X e''-[> 2](fb) 



25.6 ±2.8 ±3.4 



cJBabaAe^e- ^ J/V' ± X 2](fb) 



17.6 ±2.8 ±2.1 



Theoretical Result 




(0101(35^)10) 
3 


(0|Pi(i55-)|0) 


(0|Pi(35f)|0) 
3m^ 


as 


a(fb) 


0.243Gev3 


0.252Gev3 








as = 0.26 


f^LOiasy) = 2-26 


0.337Gev3 


0.450Gev3 








as = 0.26 


(^NLOias) = 10-92 


0.286Gev3 


0.295Gev3 


0.0321Gev3 


0.0321Gev3 


as = 0.26 


(^NLOiv^) = 3-87 


0.432Gev3 


0.573Gev3 


0.0514Gev3 


0.0514Gev3 


as = 0.26 


(^NLOias,v^) = 20.04 



16.8fb to 23.0fb when rUc varies from 1.4GeV to 1.6 GeV, as shown in Fig[3]. The calculated 
cross sections of a[e~^e~ J/ip + r^c] are listed in Table I, where aLo{as,v'^) means the cross 
section at leading order in both as and f^; <JNLO{as) means that obtained by using the short 
distance part with as corrections Q], and the long-distance matrix elements also with as 
corrections extracted from both J/ip ^ e+e^ and rjc 27 (not from J/ip ^ e~^e~ alone as 
was extracted in [9]) with all at leading order in t>^; <Jnlo{v^) means that obtained with 
corrections but without as corrections; and crNLO(asy) means that obtained by combining 
both as and v"^ corrections, where the matrix element are taken from Eq.|20j. 

We see that with both QCD radiative corrections and relativistic corrections, the dis- 
crepancy between experiment and theory for e~^e^ J/ip + rjc could be largely resolved. 
Our result is consistent with 1^. In our approach the long-distance matrix elements are ex- 
tracted from experimental data of J/ip and rjc decays, and therefore are model- independent. 
The main uncertainties may come from the higher order corrections and the errors in the 
measurements. The relativistic effects on the double charmonium production estimated in 
our approach are milder than some results obtained by using the light-cone methods. 
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Lpped gia.plis 



FIG. 4: Feynman diagrams for e"*" + e ^* ^ J /ij: + cc. 
IV. RELATIVISTIC CORRECTIONS TO e+e ^ J/i) + cc 

The long distance matrix elements in process e^e~ J /ip + cc are the same as in process 
e^e~ ^ J /ip + rjc, so we just use the result above, and only give the detailed calculation for 
the short distance part. 

A. Short distance part 

There are four Feynman diagrams in the process e~^e~ {CC)-as^ + CC, which are shown 
in Fig[4]. As in the last section the amplitude can be expanded in terms of the quark relative 
momentum and reads 

M(e+e- ^ (CC),s, + CC) = (^)"''A(q^,) = 
(f)"^(-1(0)+<|l| +■■•). (22) 

where 

^(?^) = E E4' ^' ^^1^' ^^)(3,^; 3, j|l)A(e+e- ^ + q^)C^,,{^ - q^) + CC\22,) 

A1A2 ij 

At leading order of f ^ 

|M|2 = |^A(0)A*(0) + \q^^q'^A^pA*{Q) + ^g^g^/i:^/l(0), (24) 
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where A^p = ^|^, and A^^ = Then 



2Nrm, 



■(A(0)^*(0)) 



(25a) 



,d{^A{Q)A*m 



2Nrmr 



+ g«n.^(P,/^)A*(0) + Vn.^(P,/^)A(0)) 



(25b) 



The expressions of A{0)A*{0) and n"P(AagA*(0) + A;^^(0)) can be written in terms of 
some dimensionless variables Zi — 2Ei/y/s, 'qi — 2pi/^/s, Xi — cosOi and Si — 2mi/^/s. Here 
y/s is the total energy in the center of mass frame, pg, P4, are the four-momenta of the 
final state J/ip, charm quark and anticharm quark respectively, mf = pf , and 9i is the angle 
between state i and the electron. The scalar productions are expressed as follows: 

5 5 5 

P1-P3 = ^(^3 - g3a:;3); P2.P3 = ^(^3 + g3a;3); P4-P5 = -(4 - 4^3 + ^3 - ^4 - S^); 

s s 
Pi-Pa = -^{za - {q-x- - qzX^)/2)- p2.pi = ^(^4 + {q-X- - q^x^)l2)- 

5 5 
Pl-Pb = ^(^5 + {q-X- + q3X3)/2)] P2-Pb = ^(^5 - (g-2^- + q3X2)/2)] 

Pi-Pi = ^(4 - ^3 - 4^5); P3.P5 = ^(4 - 52 - 4^4); P1.P2 = |; (26) 

where z_ = z^-z^, q^ = \ql-q^\ = ,/4 - 4^3 + zi + 5'^ - 251 - 2^5, = l^sl 
x_ = cos^_, and 9_ is the angle between ql and the electron. 



\/^3 - 5|, 



B. Numerical result 



In principle, the phase space calculation also contains relativistic corrections (e.g., by 
inclusion of the binding energy in a meson) . But for simplicity, we calculate the phase space 
integration by assuming ruj/^ = 2mc (the effect due to this simplification is very small). 
Then the three-body final state phase space expressed in terms of the variables defined above 
becomes 

5 



1=3 



rf$3 = {2Ti f 6^{pl +p2-p3-pA- p5) Yl 

s dzsdxsdz-dw 
~ 32(27r)4 ^(13 X2)(i _ ^2) 



d^Pi 
2K- 



(27) 
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where 



The hmits of those variables are 



K= '-^^-''\ (28a) 
w = x^+ Kxs. (28b) 



S<Z3<1, (29a) 
- 1 < 2:3 < 1, (29b) 



;i - K^){1 -xl)<w< ^J{l - K^){1 - xi), (29d) 
where 5 = 4m d y/s. 

The expressions of Mq and Mi are lengthy, so we only give the expression for the differen- 
tial cross section. The differential cross section of unpolarized J/ip production in association 
with cc in e~^e~ annihilation can be expressed as 

— ^^(e+e- ^ 7* ^ ^ + X) = S{E,){1 + acos^^a), (30) 
aL^a cos 03, 

where is the energy of J/ip, and 6^ is the angle between J/tp and the electron. When 
expressing the above equation by using the new parameters it becomes 

[e^e-^r^i' + X) = S{zs){l + a{zs)xl). (31) 



dz?,dx?, 



The leading order result of S'(^3) and a{z^) are in agreement with those in |17l . |21| . We 
give the next-to- leading order result of ^(za) and a (-23) in the Appendix A. Using the long- 
distance matrix elements in Eq.[TS] (without NLO corrections in the charmonium decay 
widths), the leading order cross section of + e~ J/ip + cc with = 1.5 Gev, ^/s = 10.6 
Gev, as = 0.26 is estimated to be 110 fb,*^ and the next-to-leading order f ^ correction is only 
0.42 fb, which gives about less than a half percent enhancement. As a result, the relativistic 
correction to e"*" + ^ J/ip + cc is found to be very small and negligible, in contrast to 
the exclusive double charmonium production process e+ + e~ —>■ J/ip + i]c, in which the 
relativistic correction is very significant. In fact, we find that in + e~ J/'^P + cc the 



b 



This value is smaller than 148 fb given in 21|, because here a smaller value of the wave function squared 



at the origin is extracted from the J/ip data than that used in 21| from potential model calculations. 
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FIG. 5: i?^, as a function of y/s. Here is the ratio of the correction at order to the leading 
order result for the cross section of e+ + e~ — > J/V' + cc. 

short distance part Nim^ at order v"^ is much smaller than the leading order term iVo, i.e., 
^ 1. Here Ni in Eq. [25b] is the sum of two terms, and they both are small and have 
different signs, so their sum becomes even tiny. So, the tiny effect of relativistic corrections 
on the rate of + e~ — > J/V' + cc is due to the smallness of the short-distance coefficient 
correction, regardless of the long-distance matrix elements. 

When y/s is larger, the ratio of the correction at order f ^ to the leading order contribution 
even changes sign from positive to negative values, but its value is always small. This is 
shown in Fig[5], where the values of the parameters are the same as used above. 

If we use the enhanced matrix elements in Eq. [20] , which are obtained by including NLO 
Us corrections in the charmonium decay widths, for ttlc = l.SGev, ^/s = 10.6Gev, = 0.26, 
the leading order result is about 214fb, and the relativistic correction at order v"^ is only 
0.67fb. 



C. Discussions 

Since for e"^ + J/ip + cc the conffict between experiment and theory at leading 

order in both a, and f ^ is quite serious, a number of attempts have been made to solve this 



problem. In [22(] the authors considered the production of J/'ip + cc through two photons, 
and found that the contribution of two photon process may be comparable to that of one 
photon process when y/s is large, say, larger than a/s = 20Gev, but at a/s = 10.6Gev, it is 
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only 15% of the one photon contribution. Authors in 3l| considered the process when the 



mass of the charm quark can be neglected. In 



32| the result in factorization method was 



compared with the duality method. Other suggestions can also be found in e.g. [33|, |34 |. 



Despite of these efforts, a satisfactory resolution is still needed. Recently, the authors of 



ref.j23|] have calculated the next-to-leading order QCD corrections to the direct J / ip + cc+ X 
production in e"*" + e~ annihilation at ^/s = 10.6GeV, and with rric = 1.5GeV and as = 0.26 
their result is 

(Tdirectie^e- ^ J/^ + cc) ^ 0.33pb, (32) 

It enhances the leading order result by a factor of 1.8. Further including the contributions 
from 'ip{2S) and other states, the cross section of prompt J/ip + cc production becomes 

CTpromptie^ e~ J/ijj + CC + X) ^ O.Slpb. (33) 

It is about 60% of the Belle value in Eq. [2]. This result shows that the QCD radiative 
corrections to e'^ + ^ J/tp + cc is essential. 

Here, when we further consider the relativistic corrections at order f^, we find that in 
contrast to the J/ip + rjc production, the relativistic corrections to J/ip + cc production 
are very small. The relativistic corrections can only improve the leading order result by 
enlarging the long distance matrix elements, and have little effect on the short distance 
part. As a result, the relativistic corrections to + —>■ J/ip + cc are tiny and negligible. 

The physical reason for the difference between exclusive and inclusive processes largely 
lies in the fact that in the exclusive process e~^e~ — > J/tp + rj^ the virtuality of the gluon 
that converts into a charm quark pair takes its maximum value | in the nonrelativistic 
limit (with zero relative momentum between quarks in charmonium), and introducing the 
relative momentum can substantially reduce the gluon virtuality, and hence enhance the 
short- distance coefficient and the cross section; whereas for the inclusive process e"*" + e~ — >■ 
J/ip + cc, the charm quark fragmentation c—^J/ip + clc—^J/tp + c) is significant, in which 
the quark relative momentum has little effect on the virtuality of the gluon. Therefore, our 
result for the next-to-leading f ^ corrections may indicate that the relativistic correction is 
not a good direction to solve the problem of discrepancy between theory and experiment in 
J/ip + cc production. Further studies for QCD radiative corrections and other possible new 
mechanisms related to the double charm production could be more useful. 
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V. SUMMARY 



In this paper, we have studied the relativistic corrections at order v"^ to the double-charm 
production processes e^e~ —>■ J/tp + rjc and e~^e^ —>■ J/ip + cc atB factories in the framework 
of non-relativistic quantum chromodynamics. The short-distance parts of production cross 
sections are calculated perturbatively, while the long-distance matrix elements are extracted 
from experimental data for J/tp and rj^ decays up to errors of order f^, and therefore are 
model-independent. The main uncertainties may come from the higher order corrections 
and the errors in the measurements. Our results show that the relativistic correction to the 
exclusive process e~^e~ ^ J/ip + t]c is significant, which, when combined together with the 
next-to-leading order as corrections, could resolve the large discrepancy between theory and 
experiment. It can be clearly seen from TABLE I that for the cross section the relativistic 
correction alone gives an enhancement factor of 1.7 while the combination of relativistic 
correction with QCD radiative correction results in a much larger enhancement factor of 
9. This conclusion is consistent with [l^, in which, however, the long-distance matrix 
elements are estimated by using potential model calculations. The relativistic effects on 
the J/iprjc production estimated in our approach are milder than some results obtained by 
using the light-cone methods. On the other hand, for the inclusive process e~^e~ J/ip + cc 
we find that the relativistic correction is tiny and negligible, and therefore not helpful in 
resolving the discrepancy between theory and experiment. The physical reason for the 
above difference between exclusive and inclusive processes largely lies in the fact that in 
the exclusive process the relative momentum between quarks in charmonium substantially 
reduces the virtuality of the gluon that converts into a charm quark pair, but this is not 
the case for the inclusive process, in which the charm quark fragmentation c — > J/ip + c 
(c J/ip + c) is significant, and QCD radiative corrections can be more essential. Further 
studies are needed to understand the large ratio of the cross section of e~^e~ J/ ip + cc to 
the cross section of e~^e~ — > J / ip + anything, which is measured to be in the range of 0.6-0.8 
by Belle. 
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APPENDIX A 



The formulas oi S{z3) and ^(Za) are given as follows: 



- (98304 + 98304 z^^ + 696320 - 6963200 z^^ + 14895104 z^^ - 13801472 ^3^°+ 
6104576 Z2^^ - 1221632 ^3^^ + 93596 zj}^ - 737280 ^3^ 5^ + 3047424 ^3^ 5^ - 6987776 z^^ (5V 
15003648 ^3^ 5^ - 22528000 ^3^ 5^ + 17085440 ^3^ 5^ - 4745728 ^3^ 5^ - 411904 ^3^° 5^ 
299584 ^3^^ 5^ - 26176 ^3^^ 5^ + 737280 z^ 5^ - 4030464 z^^ 5* + 7950336 z^^ 5^- 
7868416 ^3^ 5^ + 3901440 z^^ 5^ + 2526208 z^^ 5^ - 6324224 ^3^ 5^ + 3657984 ^3^ 5^- 
540096 ^3^ 5^ - 42816 ^3^° 5^ + 3632 z^^^ 5^ + 442368 5^ - 1124352 z^ 5^ + 387072 z^^ 5^+ 
948224 Z3^ 5^ - 1466368 5^ + 1291008 ^3^ 5^ - 68864 zs^ 5^ - 493440 ^3^ + 99136 ^3^ 
34552 ^3^ 5^ - 664 zs^^ 6^ + 211968 6^ - 489984 ^3 5^ + 357888 ^3^ - 14336 zs^ 5^- 
365056 ^3^ 5^ + 357696 ^3^ 5^ - 42464 z^^ 5^ - 23104 z^^ 5^ - 7356 ^3^ 5^ - 138 ^3^^ 5^+ 
39168 - 93696 z^ 5^° + 95232 ^3^ 5^° - 20224 ^3^ - 52576 zs'^ 5^° + 23392 Z3^ 5^°+ 
3904 ^36 5^° + 432 ^3^ (5^° + 45 ^3® + 3456 S^^ - 9216 ^3 + 8352 ^3^ 

384 z,' 6'' - 2520 z^' 6'' + 192 ;.3' S'' - 90 ^3' 5^^) to ^ ~ ^2 J 

3(2 - 2:3)2 y 52 - 4^3 + 4 

- 192(256 - 512 Z3 + 320 zs^ - 64 ^3^ + 16 4 - 16 ^ + 4 ^ - 128 5^ + 64 Z3 5^+ 

240 52 -224 4 52 + 20 4 5^ + 16 4(^^ + 485^-96^35^ + 2445^ + 2445^ + 345^' 

5^4(1 ~ -^3) arctan 



52 - 4^3 + 4 

(8192 ^3^ - 102400 4 + 262144 4 - 282624 4 + 146944 4° - 35584 4^ + 3072 4V 
69632 4 5^ - 299008 4 + 542720 4 - 564224 4 5^ + 331264 4 - 87808 4 
1664 4 52 + 1472 4° 52 + 336 4^ - 61440 ^3 5^ + 299008 4 - 486400 4 5^+ 
388096 4 - 180480 4 + 65280 4 - 19008 4 - 2880 4 + 3264 4 f^^- 
288 4° - 36864 5^ + 72192 Z3 5^ - 16896 4 - 64256 4 + 65664 4 f^^- 
15360 4 5^ - 6240 4 + 304 4 + 2216 4 + 6 4 - 8448 5^ + 14336 Z3 5^- 
6144 4 5^ - 4864 4 5^ + 1440 4 + 2624 4 5^ - 1056 4 + 15 4 
1152 5^° + 3072 Z3 5^° - 2592 4 5^^^ + 1344 4 5^° - 696 4 5^° + 240 4 5^" - 30 4 5^°) 

.^VT^m^^^^^'^^+V/ijtgMzjl). (Ai) 



a[zs)^[zs) - 35,354^4^/1^(2 - ^3)^(^1 - 5^) 

- (98304 zs^ + 98304 zl + 696320 - 6963200 ^3^° + 14895104 zl^- 
13801472 + 6104576 4^ - 1221632 4^ + 93696 - 835584 ^ 5^ + 3342336 ^ 
11059200 4 + 31518720 4 5^ - 49162240 4 5^ + 37194752 4° _ 11941888 4^ (^^^ 
816896 4^ + 153152 4^ 5^ - 26176 4^ - 1867776 4 + 6881280 4 f^^- 
4657152 4 S'^ - 19660800 4 + 43751424 4 - 31881216 4 + 6889984 4 5^+ 
366080 4° 5^ + 119808 zf 5^ + 31744 zf 5^ - 976 zf 5^ + 245760 ^3 5^ - 1622016 4 5^+ 
3229696 4 - 1411072 4 5^ + 2289664 4 5^ - 12958720 4 5^ + 16642816 4 f^^- 
7645696 4 5® + 1814848 4 5^ - 524672 4° 5^ + 16520 4^ - 1240 4^ + 147456 5^- 
522240 Z3 + 967680 4 - 2955776 4 + 6496768 4 - 5705984 4 5^ + 1570816 4 f^^- 
345536 4 + 177184 4 5^ + 68040 4 + 3420 zf 5^ + 294 4^ 5^ + 46080 5^° + 10752 ^3 5^°- 
254208 4 5^° + 92672 4 5^° + 144640 4 5^° + 86976 4 + 9920 4 5^° - 81632 4 5^°- 
996 4 5^° - 2070 4 (5^° + 45 4° 5^° + 16128 S^^ - 18432 ^3 S^^ - 15744 4 (^^' + 65408 4 5^'- 
116416 4 (^^^ + 61312 4 5^^ - 4136 4 + 1800 4 (^^^ + 45 4 + 3456 5^^ - 9216 ^3 5^^+ 

8352 4 + 384 4 - 2520 4 S'' + 192 4 - 90 4 5^^) ^^/^^ ~ . / " ^ 

^ 3(2 - zs f y 52 - 4^3 + 4 

+ 192(256 - 512 Z3 + 320 4 - 64 4 + 16 4-164 + 44" 128 + 64 Z3 5^ + 240 4 f^^- 
224 4 52 + 20 4 5' + 16 4 5' + 48 5^- 96 Z3 5^ + 24 4 5^ + 24 4 5^ + 3 4 5^) 

r2 _ A2 

(5S^(1 - ^3) arctan + (32768 4 - 90112 4 + 20480 4 + 180224 4°- 

251904 4^ + 140800 4^ - 35072 4^ + 3072 4^ - 258048 4 + 806912 4 5^ - 845824 4 
87040 4 5^ + 484864 4 5^ - 382208 4° + 116096 4^ - 12608 4^ - 48 4^ 
49152 4 5^ - 286720 4 + 642048 4 - 845824 4 5^ + 871680 4 - 684544 4 5^+ 
341824 4 5^ - 90112 4° 5^ + 9712 4^ 5^ - 288 4^ - 20480 ^3 + 122880 4 f^^- 
289280 4 5^ + 453120 4 6^ - 582144 4 5^ + 489856 4 5^ - 242112 4 5^ + 64864 4 5^- 
6352 4 5^ - 1336 4° - 138 4^ - 12288 5^ + 36352 ^3 5^ - 70400 4 + 136448 4 f^^- 
160640 4 + 127616 4 5^ - 62400 4 5^ + 12272 4 5^ + 1256 - 6 2'' 6^ + 15 4° 5^- 
768 5^° - 10752 5^° + 29568 4 5^° - 36480 4 5^° + 18368 4 5^° - 2144 4 5^° - 472 4 5^°- 
24 4 (^^° + 15 4 (^^° - 1152 5^^ + 3072 Z3 - 2592 4 (5^^ + 1344 4 - 696 4 (^^V 



240 4 - 30 4 ^-)^^^ln ^3^ TO^ + 2 y(^3^^)(4^ ^ _ 

zsVS'' - 4^3 + 4 - 2 V(l - zs)(zi - 5^) ) 



(A2) 
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